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A number of issues related to the flutter and postflutter of two-dimensional supersonic lifting surfaces are 
addressed. Among them there are the 1) investigation of the implications of the nonlinear unsteady aerodynamics 
and structural nonlinearities on the stable/unstable character of the limit cycle and 2) study of the implications 
of the incorporation of a control capability on both the flutter boundary and the postflutter behavior, lb this 
end, a powerful methodology based on the Lyapunov first quantity is implemented. Such a treatment of the 
problem enables one to get a better understanding of the various factors involved in the nonlinear aeroelastic 
problem, including the stable and unstable limit cycle. In addition, it constitutes a first step toward a more general 
investigation of nonlinear aeroelastic phenomena of three-dimensional lifting surfaces. 
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Nomenclature 

speed of sound 
semichord length 

linear viscous damping coefficients 
in plunging and pitching, respectively 
plunging displacement and its dimensionless 
counterpart (== h /b), respectively 
mass moment of inertia per unit span 
linear stiffness coefficients in plunging 
and pitching, respectively 
nonlinear stiffness coefficient in pitch and its 
normalized counterpart [Eq. (5)], respectively 
unsteady lift and moment per unit wing span, 
and their dimensionless counterparts, 

( Lab/mU ^ ) and (M a b 2 / /<*£/£,), respectively 

nonlinear moment control; linear 

and nonlinear control gains 

nonlinear restoring moment 

undisturbed flight Mach number, I/ooMw* 

and its normalized counterpart, Moo/(f^Xar a )* 

respectively 

airfoil mass per unit span 
pressure, air density, and speed of sound 
of the undisturbed flow, respectively 
dynamic pressure, jPocU^o 
dimensionless radius of gyration with respect 
to the elastic axis, ^/U a /mb 2 ) 
static unbalance about the elastic axis 
and its dimensionless counterpart, 

S a /mb , respectively 
time variable and its dimensionless 
counterpart, U^t/h, respectively 
free stream speed and its dimensionless 
counterpart, Uoo/ba > a , respectively 
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Superscripts 
\ f 


speed and frequency of flutter 

downwash velocity normal 

to the lifting surface 

transverse displacement 

elastic axis position measured from the 

leading edge (positive aft) and its 

dimensionless counterpart, *ea /b, 

respectively 

twist angle about the pitch axis 
aerodynamic correction factor 
tracing quantities identifying the structural, 
aerodynamic, and nonlinear control terms, 
respectively 

damping ratios in plunging and pitching, 
Ch/lmMh and c a /2I a o) a , respectively 
isentropic gas coefficient 
dimensionless mass ratio, m/Apb 1 
linear and nonlinear normalized control gains, 
f\/K a and fit K a * respectively 
uncoupled frequencies of the linearized 
aeroelastic system counterpart in plunging 
y/(K h fm), pitching y/{K a /I a ), and frequency 
ratio ( u)h/(o a ), respectively 


time derivative and its dimensionless 
counterpart, respectively 


Introduction 

H IGH-SPEED and high-performance combat aircraft perform 
aggressive maneuvers that can result in significant reductions 
in the flutter speed. Moreover, the tendency to increase structural 
flexibility and maximum operating speed increases the likelihood of 
flutter within the aircraft operational envelope. This can jeopardize 
aircraft performance and dramatically affect its survivability. To 
prevent such events from occurring, two principal issues need to be 
addressed: 1) increase of the flutter speed without weight penalties 
and 2) investigation of the possibilities of converting the unstable 
limit cycle into a stable limit cycle. The successful accomplishment 
of the second issue will permit the crossing of the flutter boundary 
without danger of a catastrophic failure. In such a case, however, 
structural fatigue becomes a concern. 

Before addressing these issues, the search for the aeroelastic in- 
stability of lifting surfaces encompasses two basic problems. One 
of these, based on the linearized aeroelastic equations, allows de- 
termination of the flutter boundary. The second one, based on the 
nonlinear approach to the aeroelastic problem, allows determination 
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of the character of the flutter boundary. In this sense, the flutter 
boundary can feature either benign or catastrophic behavior. 

Because of the necessity of avoiding flutter and/or flutter-related 
airplane performance restrictions, it appears that determination of 
both the flutter boundary and of its character, that is, catastrophic or 
benign, and the possibility of controlling both of these present con- 
siderable practical importance. The goal of the control is to expand 
the flight envelope without weight penalties by increasing the flutter 
speed and to convert the catastrophic flutter into benign flutter. 

The concept of catastrophic and benign types of flutter can be 
found in the specialized literature under different connotations that 
depend on the particular approach of the problem. The terminology 
of benign or catastrophic flutter 1-4 is synonymous with that of stable 
and unstable limit-cycle oscillation (LCO), 5-11 also referred to in 
the literature as supercritical and subcritical Hopf bifurcation 12 (also 
Refs. 7, 13, and 14), respectively. The various terminologies related 
to the character of the flutter boundary and a few sources where 
these can be found are shown in Table 1. These terminologies are 
used throughout the paper. 

In this study, the issues related to both the increase of the flutter 
speed and the character of the flutter boundary, as well as of their 
control, will be addressed. In the aeroelastic governing equations, 
the various nonlinear effects on which basis is possible to analyze 
the character of the flutter boundary will be incorporated. An ac- 
tive control methodology capable of expanding the flutter bound- 
ary and of converting the unstable LCO into a stable LCO will 
be implemented. The nonlinearities to be included in the aeroelas- 
tic model can be structural, that is, arising from the kinematical 
equations, 7-9,11 physical, that is, those involving the constitutive 
equations, 2-4,1516 or aerodynamic appearing in the unsteady aero- 
dynamic equations. 2-4,6,17 Their contribution can be beneficial (be- 
nign flutter boundary) or detrimental (catastrophic flutter boundary). 
A discussion of this issue in the context of the panel flutter may be 
found in Refs. 2-5. 

The nature of the LCO, which provides important information on 
the behavior of the aeroelastic system in the vicinity of the flutter 
boundary, can be examined by the nature of the Hopf bifurcation 12 of 
the associated nonlinear aeroelastic system. 7,13,14 Figure 1 presents 


Ihble 1 Terminologies of the dynamics of nonlinear 
aeroelastic system 


Terminologies 

Selected references 

Benign/catastrophic 

Bolotin 1 and Librescu 2-4 

flutter boundary 
Stable/unstable LCO 

Dowell, 5 Friedman and Hanin, 6 

Supercritical/subcritical 

Lee et al., 7 Lee and Kim, 8 Mei 9 
Morino, 10 Strganac et al., 11 
Holmes, 13 Lee et al., 7 

H-B 12 

Mastroddi and Morino 14 



Fig. 1 Character of the flutter boundary in the terms of LCOs ampli- 
tudes. 


several pertinent scenarios; V = V F defines the flutter boundary that 
can be determined via a linearized analysis. The nonlinear approach 
to the problem enables one to determine the aeroelastic behavior in 
the vicinity of the flutter boundary. As a result of the nonlinear anal- 
ysis, one can determine the aeroelastic behavior for a flight speed 
lower than the flutter speed V F (curve 1), that is, for V < V F , where a 
subcritical aeroelastic response is experienced. For V > Vp, the sys- 
tem can exhibit either a stable LCO (supercritical Hopf bifurcation 12 
(H-B), curve 2), or an unstable LCO (subcritical H-B, curve 3). 

In this paper, a general approach to the problem of the stability 
of the LCO of supersonic/hypersonic two-dimensional lifting sur- 
faces is addressed. This methodology enables one to accomplish a 
parametric study over a large number of parameters that character- 
ize the aeroelastic system. 17 Literature dealing with the problem of 
the determination of the flutter boundary of a supersonic/hypersonic 
wing section and on the nature of the LCO in the presence of both 
structural and aerodynamic nonlinearities is quite scarce. 2-4 

Nonlinear Model of the Wing Section 
Incorporating Active Control 

The aeroelastic governing equations of controlled wing section 
featuring plunging and twisting degrees of freedom, elastically con- 
strained by a linear translational spring and a nonlinear torsional 
spring exposed to a supersonic/hypersonic flow field are 18 

mh(t) + S a a(t) + c H h(t) + K h h(t) = L a (t) ( 1 ) 

S a h(t) + l a <x(0 + c a a(t ) + M a = M a {t) - M c ( 2 ) 

where h(t) is the plunging displacement (positive downward), a (t) 
is the pitch angle (positive nose up), and the superposed dots denote 
differentiation with respect to time t. Moreover, in Eq. (2) 

M a =K a a(t) + S s K tt a\t) (3) 

represents the overall nonlinear restoring moment that involves both 
the linear and the nonlinear stiffness coefficients, K a and K a , respec- 
tively. The tracer 8s in Eq. (3) can take the value 1 or 0 depending on 
whether the nonlinearity is included or ignored, respectively. Within 
a linear model (5s =0), M a reduces to K a ot(t ). The nonlinear coef- 
ficient K a in Eq. (3) can assume positive or negative values. Positive 
values of K a account for hard structural nonlinearities, whereas neg- 
ative values of K a account for soft structural nonlinearities. (Notice 
that this nonlinearity appears in the present case in the equation 
relating the restoring moment with the pitch angle and that it has 
the character of a constitutive equation. For this reason, it would be 
more appropriate to refer to these as physical nonlinearities.) The 
active nonlinear control can be represented in terms of the moment 
M c in a similar functional form as 

M c = fMO + Scfi^U) ( 4 ) 

In Eq. ( 4 ) f\ and / 2 are the linear and nonlinear control gains, 
respectively. Within a linear active control methodology, the tracer 
assumes the value 5c = 0. Reducing the aeroelastic equations to 
dimensionless form, we define the parameter B that represents a 
measure of the degree of the structural nonlinearity of the system 
and two normalized linear and nonlinear control gain parameters \jr x 
and i/f 2 , respectively, as 


B = K a /K a (5a) 

= fi/K* (5b) 

^2 = fi/K a (5c) 

Corresponding to B < 0 or B > 0, the structural nonlinearities are 
soft or hard, respectively, whereas for 5=0, the system is struc- 
turally linear. 

The nonlinear unsteady aerodynamic lift and moment, from 
piston theory aerodynamics (PTA), 19,20 defines pressure on the 
upper and lower faces of the lifting surface as p(x,t) = 
PoJl + v z (k - lV^oo] 2 *^* -11 , where =* J p 00 /p 00 ; v z is the 
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downwash velocity normal to the lifting surface expressed as 19 v z = 
— (dw/dt -}-f/ 00 9iu/3x)sgnz; w is the transversal displacement of 
the two-dimensional lifting surface, w(t) = h(t) + a(0(* -*ea); 
and sgn z is the sign distribution that assumes the values 1 or — 1 for 
z > 0 and z < 0, respectively. In addition, jcea = bxo is the stream- 
wise position of the pitch axis measured from the leading edge 
(positive aft). 

Retaining in the binomial expansions of p(x, t ), the terms up to 
and including ( v z /a oo) 3 , yields the pressure formula for the PTA in 
the third-order approximation 2,3 ’ 21 

p/Poo = 1 +Kyv z /a 9 o +*(* + l)/4(yr z /aoo) 2 

+ *(* + l)/12(yv-/aoo) 3 (6) 

The aerodynamic correction factor y = AfooA/(^oo — 1) enables 
one to extend the applicability of the PTA to the low supersonic 
flight-speed range. 21,22 Equation (6) is valid as long as the transfor- 
mations through compression and expansion are considered to be 
isentropic, that is, as long as the shock losses would be insignifi- 
cant (low-intensity waves). On the other hand, a more general pres- 
sure expression, obtained from the theory of oblique shock waves 
(SWT), that is valid over the entire supersonic/hypersonic range was 
obtained in Refs. 21 and 22, and it was used in aeroelastic analyses 
in Refs. 2-4. It is given by 

p/Poo = 1 4* KyVz/aoo + tc(jc + l)/4(y 

+ K(K + \) 1 /'i2{.YV.Ja ca ? (7) 

With the exception of the cubic terms, Eqs. (6) and (7) resemble 
each other. This is explained by the entropy variation appearing 
in the pressure expansion, beginning with the third-order terms. In 
contrast to Eq. (6), Eq. (7) encompasses additional features in the 
sense of 1) taking into account shock losses that occur in the case 
of strong waves, 2) being applicable over a wider range of angles of 
attack (a < 30 deg) and Mach numbers (Af > 1 .3) (Refs. 2 1 and 22), 
and 3) being applicable to Newtonian speeds ( M — ► oo, y — > 1). 
Comparison of results showing the unstable and stable LCOs using 
the PTA and SWT will be presented next. 

The two coefficients of the cubic terms in the two equations differ 
by 10% for k = 1.4, and so, for a more accurate prediction of the 
character of the flutter instability boundary, it should be included 
(Ref. 20). However, within a linear stability analysis, the flutter 
speed evaluated via these two expressions by SWT and PTA does 
not exhibit any differences. 

A comparison of the flutter speed vs flight Mach number obtained 
from the PTA and SWT, including and discarding the correction 
factor y, is shown in Fig. 2. In addition, in the same figure, the 
flutter boundary obtained via the use of the linearized supersonic 
unsteady aerodynamics as provided by Garrick and Rubinow, 23 is 
also supplied. In the low -supersonic flight- speed regime, the PTA 
and SWT with the corrective term provide a rather good agreement 
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Fig. 2 Comparison of the predictions of the flutter speed vs the flight 
Mach number when using PTA, SWT, and the exact unsteady super- 
sonic aerodynamics: p = 100, Xa = 0*25? <*> = 1*2, r a - 0.5, ^ a ~Ck= 0, 
and *0 = 0.5. 


with the flutter predictions reached via the supersonic flow theory, 23 
and as a result, this correction should be included. At the same time, 
for higher supersonic Mach numbers, the differences in the flutter 
predictions based on the indicated aerodynamic theories practically 
disappear. In the next developments, unless otherwise stated, PTA 
will be applied. A comparison of the predictions of the benign and 
catastrophic character of the flutter boundary, based on these two 
aerodynamic theories, will be shown subsequently. 

When the case of the flow on both surfaces of the airfoil with the 
speed = U Q 0 is considered, from Eq. (6) the aerodynamic 

pressure difference Sp can be expressed as 

&p\n\ = (4^/Af 0C ,)y[(tti, / /i7 (K + w, x ) 

+ (l+«c)/l2 y 2 M 2 m (w„/U 00 + w, x ) 3 ] (8) 

In the next developments, the nonlinear aerodynamic damping in 
Eq. (8), that is, the terms associated with (w,,) 2 and (tu,,) 3 , will be 
discarded, and consequently, the cubic nonlinear aerodynamic term 
reduces to the (w,* ) 3 only. The study of the implications of the non- 
linear aerodynamic damping on the nature of the LCO constitutes 
an important problem, which is not addressed in this paper. 

Next, the nonlinear unsteady aerodynamic lift L a (t) and moment 
M a (t ) per unit wing span can be obtained from the integration of the 
pressure difference on the upper and lower surfaces of the airfoil: 

fib 

L a (t)= / 8pdx (9a) 

Jo 

f 2 * 

M a (t) = - I 8p(x — .x'ea) dr (9b) 

Jo 

The final expressions can be cast in compact form as 

Lad) = — (M/ooPoo/3M 00 )}' { XlUao^d) 

+ 5,M^i/ 0O (l +K)j/V(f) + 12 [h(t) + (b - jrEA)a(f)]} 

(10a) 

Mad) = (6(/ooPoo/3Afoo)y { HU^ib - *ea )<*d) 

+ S A Ml 0 U 00 (b - ^ea)( 1 + K)vV(r) + 4[3 (b - XEA )h(t) 

+ (4b 2 — bbx VA + 3.r I 2 - A )a , (r)] } (10b) 

where <5,4 is a tracer that is set equal to 1 if the aerodynamic nonlin- 
earity is included or set equal to 0 if the aerodynamic nonlinearity 
is ignored. 

As a result, the governing equations ( 1 ) and (2) considered in con- 
junction with Eqs. (10a) and (10b) feature inertial and aerodynamic 
coupling. Using the dimensionless time r = U^t/b, the system of 
governing equations can be expressed as 

$'(*) + XaCt"(T)+2 C*(d>/V)f (T) + (&/V) 2 $(t) = /„( T) (11) 

(x«/ r «H"< r > + a "< r > + (2f./V)a'(r) + 1 /V 2 a(r) 

+ 1/ V 2 Bot\r) =m a (T) - ^/v 2 a( r) - f 2 /V 2 a\z) (12) 

In these equations, the dimensionless aerodynamic lift and moment 
are represented as 

l a (r) = -{Y/n^M x ){\2a(z) + & A M 2 x (.\+K)Y 2 a\z) 

+ 12K'(T) + (b-Jt EA )/ba , (T)]} (13a) 

m a (z) = ( Y /\2hM x )(\ ft*) { 12(6 - XEA )/ba(,z) 

+ S A Ml 0 (b - -*ea)/6(1 + k)y 2 o 2 (, r) + 4[3(6 - ^ea)/6^'(t) 

+ (46 2 - 6 b XEA + 3.4,) j b 2 a'(x )] } (13b) 

where £ = h/b is the dimensionless plunging displacement and the 
primes denote differentiation with respect to dimensionless time r. 
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When the procedure developed by Bautin 24 and Lyapunov 25 is 
used, pertinent conditions defining the character of the flutter bound- 
ary (benign or catastrophic), can be determined. These conditions 
are expressed in terms of the sign of the Lyapunov first quantity 4 
L(V f ) determined on the flutter boundary. 2-4 24 Specifically, the 
inequalities L(V F ) < 0 and L(V F ) > 0 define the benign (supercrit- 
ical) and catastrophic (subcritic al) nature of the flutter boundary. 
The application of Bautin’s procedure 24 requires that the character- 
istic equations obtained on the flutter boundary exhibit either one 
root or two roots that are purely imaginary. These conditions are 
equivalent to the H-B theorem. 12 

The Lyapunov first quantity 4 L(V F ) corresponding to the nonlin- 
ear flutter of the wing section in a supersonic/hypersonic flowfield is 
derived next and is used to determine the conditions that characterize 
the nature of the flutter boundary. 

The governing equations (11) and ( 12) are converted to a system 
of four differential equations in the form 2-4,24 

,4 

— - = + Pj(x U X 2 ,Xs 1 X4), j - 1,4 (14) 

dr 

m = \ 

The functions Pj(x\ , x 2 , * 3 , * 4 ) include the structural, aerodynamic, 
and nonlinear control terms that can be cast as 

Pj(x u x 2 ,x i , x A ) = ^ afxf + 2 ^ ajfxiX, + a$ xf 

i = 1 i,/«=l 1 = 1 

<M0 


+ 3^xfx l + 6 ^2 a\ilx iX ,x k 


For the present case, Eqs. ( 14) and (15) reduce to a state-space form: 

4*1 , 1 *.. 


^ = aj 3 , *i(r) + a 2 > x 1 (z) + S A a 222 x 2 (z) + & s a 2 ^xl{z) 
dr 

+ Sca^xlir) + a 2 } x } (z) + a®x 4 (r) 

^ = ai 4 ) x,(r) + a 2 ) x 1 (z) + S A a'^ 2 xl(z) + (r) 

dr 

+ 5cO^Jc|(t) 4- a^ 4> x 3 (r) 4- a‘ 4i x 4 (r) 


where 

|=jc t , a — x 2 , £'=* 3 , a f = x 4 (16e) 

The linear active control appears in the coefficients a ^ and af\ 
whereas nonlinear control terms are included in the terms accom- 
panying the tracer <5 C . The coefficients of Eqs. (16) are given in 
Appendix A. When the streamwise position of the pitch axis co- 
incides with that of the midchord, Xq = 1 , the nonlinear aerody- 
namic terms become negligible as the aerodynamic pitching moment 
vanishes. 

When the solution of Eqs. (16) is considered in the form 
Xj = Aj exp(wt), the characteristic equation corresponding to the 
linearized system obtained on the flutter boundary is 

fc> 4 + pa? -f qa) 1 + ra> + 5 = 0 (17) 


where 

P = {Vy[4+ 3 j^ + 6.*o(Xa - 1) - 6x„] + 3 r 2 (Vy 

+ 2 Moo/xCa + 2 M x nw^)}/ (r 2 - x„) (18a) 

q — {3Af 0o /i^|_A/ ( x)M(l + <*>' + lA'i) + 2t; a {V y 4- 2M 00 /xd>£; I )J 
+ Vy\V(y + - iM^fiXa) 

+ 8 + 6M x H(bx^ h 

- 3 M^x 0 (V + 4«f*)] }/3 M> 2 K 2 (r„ 2 - xt) (18b) 


5 = arfMoo/xO + \!r x )rl + V 2 y(x 0 - 1)]/ M x pV*(r 2 - 

(18c) 

r = 1 3r 2 \2M x fiw 2 ^ a + ( 1 + fi)(Vy +2M x fiw( h )\ 

+ V yib^ioixl — 6x 0 (<o + V( h ) 

+ 2(2a> + 3Vf*)]}/3M 0O / 4 K 3 (r o J - x* 2 ) (18d) 

Three of these parameters, namely, q , 5, and r, include the linear 
control gain • This has important implications on the possibility 
of controlling the nature of the LCO and the flutter boundary. For 
steady- state motion, the equilibrium is stable if the real parts of 
all of the roots of the characteristic equation are negative. 25 26 The 
analysis of the roots of this equation leads to the Routh-Hurwitz 
(R-H) conditions which define the parameter bound for the stability 
of the system. The R-H conditions reduce to the inequalities p > 0, 
q > 0, r > 0, s > 0, and 31 = pqr — sp 2 — r 2 > 0. For the aeroelastic 
stability problems, for which the condition A| =sp/r + /r/4 > 0 
should be satisfied, the roots of the characteristic equation on the 
critical flutter boundary 31 — 0, are given by 

a)\ 2 = ±/c, o >3 4 = — e ± in where i = \/—T (19a) 

c 2 —r/p , e = +p/2 

n 2 = sp/r — p 1 / 4, n > 0 (19b) 

Equation (19a) is a necessary requirement for the determination of 
L(V f ) and constitutes a statement of the applicability of the H-B 
theorem as well. 12 * 24 

For sufficiently small values of the speed V, all of the roots of 
the characteristic equation are in the left-hand side of the complex 
plane, and the zero solution of the system is asymptotically stable. 
For the value V = V F , two roots of the characteristic equations are 
purely imaginary, and the remaining two are complex conjugate and 
remain in the left-hand side of the complex plane (critical flutter 
velocity). 

To identify the benign and catastrophic character of the flutter 
boundary, it is necessary to solve the stability problem for the system 
of equations in state-space form in the critical case of a pair of pure 
imaginary roots. The expression of the flutter boundary 31 for the 
wing section is 

SR = { W^at^yixo - 1) - y 

+ 2M 0O fit; a + 2M 00 m<w(Ta) + V y [4 + 3xjj + 6* 0 (X« - 0 

- 6 X „]} 2 - 'iMl 0 n 2 [ir 1 a [2M 00 nw 1 i; a + (1 +fi)(Vy 
+ 2M 00 ^<u&)] + V ya>[4a> + 3d>Xo + 6V t; h - 6x 0 (a) 

+ v &)]} 2 ( r a - xi) + \irl\2M x ti<b 2 t; a + (1 + V'lXV'y 
+ 2M 00 fio}^ h )\ + V ya»[4a> + 3 (bxl + 6V ^ 

-6jt 0 (*>+ y?A)]}{3r 2 (Vy +2M x tit; a +2 M x fzcb^) 

+ Vy\A + 2xl + 6jc 0 (x„ - 1) - 6 Xa ] } {3W 00 ai/- 2 [M 00 ^(1 
+ a) 2 + ^ 1 ) + 2f„(Vy + 2Af 00 /x<u?)] + Vy[SM x fj.a>^i, 

+ 6M ooittbxlSh — IMaafiXoiV + 4 <b^ h ) + V (y + SMooP 

- 3AOZX,,)]} }/27M 4 W(r; - X 2 ) 3 = 0 (20) 

Notice that this expression is general and includes the relationship 
between the flutter speed and the flutter frequency parameters eval- 
uated on 31 = 0 in terms of the basic geometrical and flight param- 
eters. In the particular case in which the structural damping ratios 
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Fig. 3 Flutter speed vs flight Mach number; influence of the linear 
control gain ipi. 


in plunging and pitching are ignored and y =>• 1 , the expressions of 
flutter frequency xf and flutter speed V F reduce to 

Xf = (co a /o>) 2 F = [r 2 + (1 — *o ) 2 + 5 

-2xa(l - -*o)]/{(l + fi)rl +<i> 2 [(l ~* 0 ) 2 +5]} (21) 

and 

y _ Uf 

F bo) a Jxf 


The dimensionless structural, aerodynamic, and control nonlinear- 
ities are defined as 


^ B Vj{rl-Xa) 

(26a) 

A/oo(l A- ac)v 3 

^NL ~ / „ v 

12 fi(r% - Xa) 

(26b) 

Cm. = ti-rr? — 7 

V 2 M - X.a) 

(26c) 


where V f is the flutter speed. Upon defining and 

having in view that, according to Eqs. (19b) c > 0, the Lyapunov 
first quantity 4 reduces to 

L(Vf) = +3»4)“21 + +32*4)4 

+ (^23^222 + -^24< 1 za) oi 22 a l2 + (^1 3^222 ^l 4a S) ai2Q, 22 

(27) 

When the expressions of Snl, A ml, and C NL [Eqs. (26)] are used 
and the procedure devised in Refs. 2-4 is applied, the benign or 
catastrophic character of the flutter boundary defined as L(V F ) < 0 
or L(Vf) > 0, respectively. These conditions can be restated as 4 

V F < V r (28a) 


xi - {^xf - i K 2 (xf - 1) 


f 1 + V' i)Xa + (1 - (1 + Vfr))] - 5 U + fi) 

(22) 


respectively. 

Equations (21) and (22) represent the dimensionless flutter fre- 
quency and flutter speed of the actively controlled system. These 
constitute the extension of the equations obtained by Ashley and 
Zaitarian , 19 to include active control. 

In Fig. 3, the dependence of the dimensionless flutter speed as 
a function of the Mach number for selected values of the feedback 
gain tyx is presented. It clearly appears that, with the increase of the 
linear control gain, an increase of the flutter speed is experienced. 
Moreover, as values of are increased, the efficiency of the active 
control is increased. 

The expression of the Lyapunov first quantity 4 is given in 
closed form in Refs. 2 and 4. For the present case, the Lyapunov 
first quantity 4 is expressed in terms of the coefficients A^] as 
(Appendix B) 

L(V f ) = (3jr/4c)(A® + A® + A® + A®) (23) 

The terms in the parentheses of Eq. (23) are expressed via the coef- 
ficients a^l appearing in Eqs. (15) as 


or 


Vf > K 


(28b) 


where 


V 2 = A,/A 2 (28c) 

In Eq. (28) the parameter A\ includes the structural nonlineari- 
ties and the nonlinear control gain parameter, whereas A 2 includes 
the aerodynamic nonlinearities. Their expressions are provided in 
Appendix C. 

In the absence of the nonlinear control and for B <0 (soft struc- 
tural nonlinearities), V r is negative. The relation V F > V r corre- 
sponds to the catastrophic flutter boundary (unstable LCO) and 
occurs for any supersonic flight Mach number. For this case, an 
unstable LCO is experienced even in the presence of the linear con- 
trol. On the other hand, for B > 0 (hard structural nonlinearities), 
the transition from benign to catastrophic flutter (from stable to the 
unstable LCO) occurs at an increased flight Mach number in the 
presence of the linear control. 

As a special case, for xq = 1, Eq. (23) reduces to the following 
form: 

L(Vp) — (5$Snl + £cCm-)|_<3l3Xa — 3 l 4 )tf 2 I + (^23 Xa ~ ^24)«22 


A ktl = ( 1 / A 0 ) [a j3 a%a 2l _a v ct^ + d y4 a^a2*«a«2 J ) (24) + (3 2 3X« - 324>a22«?2 + ( 3 i3Xa - 3 u)«i2«hJ 


whereas the coefficients of the system of Eq. (24) are 
a$ 2 = (8 s B+8cf 2 )- 2( r f_“ Xa) 

Mco(l +*)y 3 [(.*o - l)Xa +r 2 ] 

*2/x(rJ - Xa) 

= (SsSfa. + ScCmJXa — ‘5a'4nl[(.* , o — l)x<* + r «] (25a) 

_ M«>(1 + *)y 3 ( xa - l + x a ) 

12 H{rl-xa) 

= — (5j5nl + <5cCnl) + Anl(-*o — 1 + Xa) (25b) 


- 5 aAnl[( 3 i 3 '-J - 3 | 4 Xa )<*21 + {^ r l “ S 24 X «)«22 

+ ($23^ ~ %24Xa)ot22<*n + ( 2 13^ ~ SuXa)®!^] (29) 

In this case, a decrease of the influence of the aerodynamic nonlin- 
earities on the aeroelastic system is experienced. 

Stability in the Presence of Active Control 

To enhance understanding of the effect of the active control on the 
character of the flutter boundary, some explanations are provided in 
Figs. 4a and 4b. In Fig. 4a, the intersection point between the two 
curves V F and V r separates the benign flutter boundary (stable LCO) 
characterized by V F < V r from the catastrophic flutter boundary 
(unstable LCO) defined by V F > V r . For these cases; a change in 
the sign of the Lyapunov first quantity 4 L(V F ) occurs (Fig. 4b). In 
this context, the following four possible scenarios are distinguished: 
1) for V < V F , as time unfolds, a decay of the motion amplitude 
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(Stable LCO) (Unstable LCO) 


a) In the (V — Afflict) plane 



b) In the {L — Af fllght ) plane 

Fig, 4 Generic representation of the flutter boundary, L< 0 corre- 
sponds to V r > Vp and vice versa. 



0.4 0.8 U 1.6 2.0 2.4 


Fig. 5 Influence of the structural and aerodynamical nonlinearities on 
the Lyapunov first quantity 4 L for the uncontrolled wing section. 

is experienced, indicative of subcritical response; 2) for V = VV, 
the center limit cycle occurs, indicative of a periodic orbit; 3) for 
V r > V F y the LCO is stable; and 4) for V F > V r , the LCO is unstable. 

The parameters for the simulations, are chosen as /r = 100, 
=0.25, a>= 1.2, r a =0.5, & =0, * 0 = 0.5, y = h*c = 1.4, 

6,4 = <5 5 = Sc = 1 , and B — 50. 

The effect of structural nonlinearities on the character of the flutter 
boundary is studied in terms of the nonlinear parameter B [Eqs. (3) 
and (5 a)]. For the present simulation, the aeroelastic system appears 
to be characterized by a catastrophic flutter boundary in the upper 
half-plane (unstable LCO) and by a benign flutter boundary in the 
lower half-plane (stable LCO). In Fig. 5, the Lyapunov first quantity 4 
L(Vf) for cases in which the soft and hard structural and aerody- 
namic nonlinearities are involved is presented. It appears that in 
the presence of only the aerodynamic nonlinearities, Lyapunov first 
quantity becomes positive for any flight Mach number. This result 
reveals that aerodynamic nonlinearity induces a catastrophic flutter 
boundary, implying that a subcritical H-B occurs. On the other hand, 
in the presence of hard structural nonlinearities only, the opposite 
situation is experienced. At relatively moderate supersonic flight 
Mach numbers, a benign flutter boundary is encountered, which be- 
comes catastrophic as the Mach number is increased. This implies 
that for higher Mach numbers the effects of the aerodynamic non- 
linearities become prevalent. It is also shown that in the case of high 


Mach numbers, the neglect of nonlinear aerodynamic terms yields 
inadvertent predictions related to the character of the flutter bound- 
ary. Moreover, when the aerodynamic nonlinearities are discarded 
( 5 4 z= 0) for any flight Mach number, the flutter boundary is benign 
or catastrophic, depending on whether hard (B > 0) or soft ( B < 0) 
structural nonlinearities are present, respectively. 

The influence of the hard structural nonlinearities, in conjunction 
with the aerodynamic nonlinearities ( 6,4 = 1 ), for the controlled/ 
uncontrolled system is presented in Fig. 6 . The dotted lines iden- 
tify the cases in which hard structural nonlinearity are included 
(B = 50), whereas the solid lines identify the cases in which the 
structural nonlinearities are ignored ( B =0). The control acts in 
both situations toward the stabilization of the system. Also, the un- 
stable LCO that occurs when only the aerodynamic nonlinearities 
are considered can be converted to a stable LCO. 

Figures 7 and 8 show that soft structural nonlinearities (B = — 10) 
result in a catastrophic flutter boundary and that via a combined 



Fig. 6 Influence of the structural and aerodynamical nonlinearities on 
the Lyapunov first quantity 4 L In the presence/absence of linear and 
nonlinear active control gains. Aerodynamic nonlinearities retained. 



Fig. 7 Conversion of the unstable to stable LCO via nonlinear active 
control, for system encompassing structural soft nonlinearity. 



Fig. 8 Unstable LCO for systems encompassing structural soft nonlin- 
earity, aerodynamical nonlinearities, and linear active control. 






968 


MARZOCCA. UBRESCU, AND SILVA 


active control (t/^ =0.5, ^2 = lOO^i) the unstable LCO can be- 
come stable (Fig. 7). Moreover, it clearly appears that, when soft 
structural ( B < 0) and aerodynamic nonlinearities are present, the 
linear active control (Vi^i > 0, 1/^2 = 0) cannot change the character 
of the flutter boundary (Fig. 8). 

From Eq. (23), which defines the Lyapunov first quantity, 4 the be- 
nign flutter boundary is expressed in closed form. When Eqs. (28a) 
and (28b) are used, the character of the flutter boundary is exam- 
ined and has been plotted in Figs. 9-11. Each of Figs. 9-11, dis- 
play in the plane (V, Afljgh,) the benign and catastrophic charac- 
ters of the flutter boundary for the actively controlled wing section 
where X = The corresponding Lyapunov first quan- 

tity is shown in Figs. 12-14 in the plane (L, A flight ). In Figs. 9- 
15, the aerodynamic and hard structural nonlinearities have been 
included, 8& = 1 and 8s = 1; B =50. In Figs. 9-15, the transition 
from catastrophic to benign flutter is shown. Figures 9-15 help one 



Fig. 9 Stable and unstable LCOs in the presence of linear control. 



Fig. 12 Influence of the linear control on the Lyapunov first 
quantity 4 L. 



Fig. 13 Influence of the nonlinear control on the Lyapunov first 
quantity 4 L 




Fig. 14 Influence of the linear and nonlinear controls on the Lyapunov 
first quantity 4 L. 



Fig. 11 Stable and unstable LCOs in the presence of linear and non- 
linear controls. 



(U DJ L2 \A 2.9 2A 


Fig. 15 Stable and unstable LCOs with and without control in presence 
of both nonlinearities; comparison between SWT and PTA prediction. 





MARZOCCA. LIBRESCU, AND SILVA 


969 


Table 2 Stable and unstable LCOs for PTA and SWT 


10 10 L(V», 10 10 L(V», 



Aflutter 

k r , PTA 

X r , SWT 

PTA 


SWT 


1.200 

2.182 

2.263 

2.386 

-6.679 



-15.437 


1.208 

2.189 

2.256 

2.378 

-5.402 



-14.107 


1.216 

2.196 

2.248 

2.370 

-4.164 

§ 

O 

-12.817 


1.224 

2.203 

2.241 

2.362 

-2.965 

« 

a 

-11.566 


1.232 

2.211 

2.234 

2.355 

-1.803 



-10.352 


1.240 

2.218 

2.226 

2.347 

-0.676 


r 

-9.175 

-2 

1.248 

2.225 

2.219 

2.339 

0.416 



-8.033 


1.256 

2.232 

2.212 

2.332 

1.474 



-6.925 


1.264 

2.239 

2.205 

2.325 

2.501 



-5.849 


1.272 

2.246 

2.198 

2.317 

3.4% 



-4.806 


1.280 

2.253 

2.191 

2.310 

4.461 



-3.793 


1.288 

2.260 

2.185 

2.303 

5.396 

s 


-2.811 


1.296 

2.267 

2.178 

2.296 

6.304 


O 

u 

-1.857 


1.304 

2.274 

2.171 

2.289 

7.183 

i 


-0.930 


1.312 

2.281 

2.164 

2.282 

8.037 



-0.031 


1.320 

2.288 

2.158 

2.275 

8.865 



0.842 


1.328 

2.295 

2.151 

2.268 

9.668 



1.689 

JU 

1.336 

2.302 

2.145 

2.261 

10.447 



2.512 

X) 

B 

1.344 

2.309 

2.139 

2.254 

11.202 



3.312 

a 

1.352 

2.316 

2.132 

2.248 

11.935 



4.088 

0 

1.360 

2.322 

2.126 

2.241 

12.646 



4.842 



to understand the behavior of the aeroelastic system in the presence 
of linear and nonlinear active control. 

The stable and unstable characters of the LCOs via the use of the 
two different aerodynamic theories (PTA and SWT) are presented 
in Fig. 15. From Fig. 15 and Table 2, one can infer that the transition 
from the stable LCO to the unstable LCO occurs at slightly lower 
Mach numbers (less then 3%) for the PTA as compared to those pre- 
dicted by the SWT. This implies that the PTA provides conservative 
results as compared to the SWT. 

Conclusions 

An original LCO analysis was presented. In contradistinction with 
the ones generally used in the literature, where the character of the 
flutter boundary is determined from the path variation amplitude of 
displacement quantities, in the present approach this information is 
obtained via the Lyapunov first quantity. 4 

It was shown that in some circumstances, the aerodynamic and 
hard structural nonlinearities contribute in different ways to the de- 
termination of stable or unstable LCOs. It was shown that the hard 
structural nonlinearities result in a stable LCO, whereas the soft 
structural nonlinearities result in an unstable LCO. At high flight 
Mach numbers, the aerodynamic nonlinearities are primary contrib- 
utors to the destabilization of the aeroelastic system. This implies 
that with increasing the hypersonic flight speed, when the aero- 
dynamic nonlinearities become prevalent, an unstable LCO occurs 
irrespective of the presence of hard structural nonlinearities. On the 
other hand, soft structural nonlinearities (B < 0) contribute in the 
same sense, as the aerodynamic nonlinearities, to the unstable LCO. 
It is also shown that active control can be used to increase the flutter 
speed and to convert the catastrophic flutter boundary into a benign 
flutter boundary and/or to shift the transition between these two 
states toward higher flight Mach numbers. 

The issue of generating the active control moment was not ad- 
dressed. It is the authors’ belief that this can be produced via a device 
operating similarly to a spring, whose linear and nonlinear charac- 
teristics can be controlled, but additional analysis are required to 
confirm this observation. 


Appendix A: Coefficients Occurring in Eqs. (16) 

The coefficients of the aeroelastic governing system represented 
in state-space form [Eqs. (16)]: 


a 


(3) 

1 



r l - x a 2 ) 


(Al) 


( 3 ) _ Moo + ft) - V 2 v\ C *0 - !)x» + r 2 ] 

V 2 M x n(r* - Xa) 


(A2) 


$A a 222 “ 


0) _ Mood 4- k ) y 3 [(*0 - l)x« + r a ] 

12 /i.i 

( 3 ) _ „ 


I2n(r 2 - x 2 ) 


Ssa'V, = B 


V 2 {rl - xi) 


<$€<*222 ^ ft 


X« r a 


V2 ( r a ~ Xa ) 


(A3) 

(A4) 

(A5) 


<P - 


V k[(3.Xq - 6.*o + 4)x« + 3(x 0 - l)r 2 ] + 6M 0o M<«Xa''a 
3VM a cft(r 2 -xj) 


«! 4, = 


0> 2 Xa 


' V 2 {rl-xt) 


(4 ) _ v^yixo - l + x«) - MqqM^O + ’Ai) 

V 2 M xt i(r~ - x 2 ) 


Ssa 222 = —B 


V 2 (r 2 - x„) 


& ca ( ^2 = ~ft 


~ X 2 ) 


(A6) 

(A7) 

(A8) 

(A9) 

(A 10) 

(All) 

(A12) 


„(4) 


„ (4> y 3 Moo(i + y)(*o - i + Xa ) 

S * a 222 = -~r- 2 X 

12 - Xa ) 

( 4 ) V^X(Ao 1 "b Xa) ~b 
° 3 ~~ VMooM(r a 2 - X a 2 ) 

VX[(3j:q - 6*o + 4) 4- 3 (.\ 0 - l)x a ] + 6M x fir 2 { a 

3VMoo/i(^-x 2 ) 

(A13) 


Appendix B: Coefficients Intervening in Eq. (23) 
The various elements of Eq. (24), evaluated on the instability 
boundary, are expressed as 

^in = ( 1 /A 0 ) (£13^222 + a 14^222) “21 (Bl) 

A 222 ~ (1/A 0 )(a23«222 + <*24 <*222 (®2) 

A | 2 j2 = (l/Ao)(a , 23 fl 222 + ^24^222)^22^2 (®3) 

= ( 1 /A 0 )(ai 3^222 + “14^222)^12^22 ( R 4 ) 

where the adjoint of a matrix is defined as A 0 [Qf 7/ ,] _1 and A 0 = ||a,y || 
is the determinant of the matrix of coefficient a ,y. In addition, the 
parameters a tJ are 


“l2 =c 


“21 = 


“22 -C 


■p 

a? 

-p 

a? 

-p 

a? 
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+ c\ 
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Appendix C: Coefficients A j and A 2 
Appearing in Eq. (28c) 


(S s B-\-8 c ^2)r^ { 3 ^ , 3 

M = 7- v |(^I3X« - ^14)0f 2 , + ($23 Xa ~ 

if a ~ Xa) 

+ P23X« - 324)a22<*r 2 + (%UXa ~ 14>0f i 2 0f| 2 ] (Cl) 


SaM^X + *:)y 3 
\2(Xirl - Xa) 


~ Sl4(*0 ~ 1 + Xa 


- %2 a(Xq “ 1 + Xa 


- %2A(X 0 * 1 + Xa 


{{S I3 [(A 0 - l)Xa +>■«] 

)}« 2 I + {^23[(-Vo - 1)X« + r \ ] 
)} a 22 + {^23 [C*0 “ DXa + ^] 

)}o' 22 a]'2 + {^b[(^ 0 — l)x« + r a] 


-Su(x 0 - 1 +Xa)}<*i2<4-} (C2) 
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